Mon. Not. R. Astron. Soc. Q00,[T|{20|(2012) Printed 19 October 2012 (MN MeX style file v2.2) 



Early-stage young stellar objects in the Small Magellanic Cloud 

J.M. Oliveira^*, J.Th. van Loon\ G.C. Sloan^, M. Sewilo^ K.E. Kraemer^ P.R. Wood^ 
R. Indebetouw*^'^, M.D. Filipovic^, E.J. Crawford^, G.F. Wong^, J.L. Hora^, M. Meixner 
T.R Robitailleii, B. Shiao^o, J.D. Simoni^ 

^ School of Physical & Geographical Sciences, Lennard-Jones Laboratories, Keele University, Staffordshire ST5 5BG, UK 
^ Department of Astronomy, Cornell University, Ithaca, NY 14853, USA 

^ Department of Physics & Astronomy, The Johns Hopkins University, 3400 N. Charles Street, Baltimore, MD 21218, USA 

^ Boston College, Institute for Scientific Research, 140 Commonwealth Avenue, Chestnut Hill, MA 02467, USA 

^ Research School of Astronomy and Astrophysics, Australian National University, Cotter Road, Weston Creek, ACT 2611, Australia 

^ Department of Astronomy, University of Virginia, PO Box 3818, Charlottesville, VA 22903, USA 

^ National Radio Astronomical Observatory, Charlottesville, VA 22904, USA 

^ University of Western Sydney, Locked Bag 1797, Penrith South DC, NSW 1797, Australia 

^ Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, MS-65 Cambridge, MA 02138-1516, USA 

Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 

Max-Planck-Institut fur Astronomic, Konigstuhl 17, D-691 17 Heidelberg, Germany 
^'^ Observatories of the Carnegie Institution of Washington, 813 Santa Barbara St., Pasadena, CA 91101, USA 



Accepted 2012 October 18. Received 2012 October 18; in original form 2012 September 18 



ABSTRACT 

We present new observations of 34 Young Stellar Object (YSO) candidates in the Small 
Magellanic Cloud (SMC). The photometric selection required sources to be bright at 24 and 
70 /im (to exclude evolved stars and galaxies). The anchor of the analysis is a set of Spitzer- 
IRS spectra, supplemented by groundbased 3—5 jim spectra, Spitzer IRAC and MIPS pho- 
tometry, near-IR imaging and photometry, optical spectroscopy and radio data. The sources' 
spectral energy distributions (SEDs) and spectral indices are consistent with embedded YSOs; 
prominent silicate absorption is observed in the spectra of at least ten sources, silicate emis- 
sion is observed towards four sources. Poly cyclic Aromatic Hydrocarbon (PAH) emission is 
detected towards all but two sources. Based on band ratios (in particular the strength of the 
11.3-/im and the weakness of the 8.6-/im bands) PAH emission towards SMC YSOs is dom- 
inated by predominantly small neutral grains. Ice absorption is observed towards fourteen 
sources in the SMC. The comparison of H2O and CO2 ice column densities for SMC, Large 
Magellanic Cloud (LMC) and Galactic samples suggests that there is a significant H2O col- 
umn density threshol d for the detection of CO2 ice. This supports the scenario proposed by 
Oliveira et al . | ( |20 1 T ) , where the reduced shielding in metal-poor environments depletes the 
H2O column density in the outer regions of the YSO envelopes. No CO ice is detected towards 
the SMC sources. Emission due to pure-rotational — transitions of molecular hydrogen is 
detected towards the majority of SMC sources, allowing us to estimate rotational temperatures 
and H2 column densities. All but one source are spectroscopically confirmed as SMC YSOs. 
Based on the presence of ice absorption, silicate emission or absorption, and PAH emission, 
the sources are classified and placed in an evolutionary sequence. Of the 33 YSOs identified 
in the SMC, 30 sources populate different stages of massive stellar evolution. The presence 
of ice- and/or silicate-absorption features indicates sources in the early embedded stages; as a 
source evolves, a compact H ii region starts to emerge, and at the later stages the source's IR 
spectrum is completely dominated by PAH and fine-structure emission. The remaining three 
sources are classified as intermediate-mass YSOs with a thick dusty disc and a tenuous enve- 
lope still present. We propose one of the SMC sources is a D-type symbiotic system, based on 
the presence of Raman, H and He emission lines in the optical spectrum, and silicate emission 
in the IRS-spectrum. This would be the first dust-rich symbiotic system identified in the SMC. 

Key words: astrochemistry - circumstellar matter - galaxies: individual (SMC) - Magellanic 
Clouds - stars: formation - stars: protostars. 
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1 INTRODUCTION 

Star formation is a complex interplay of various chemo-physical 
processes. During the onset of gravitational collapse of a dense 
cloud, dense cores can only develop if heat can be dissipated. The 
most efficient cooling mechanisms are via radiation through fine- 
structure lines of carbon and oxygen, and rotational transitions in 
abundant molecules such as water (e.g., |van Dishoeck|2 004V These 
cooling agents all contain at least one metallic atom. Furthermore, 
dust grains are crucial in driving molecular cloud chemistry, as 
dust opacity shields cores from radiation, and grain surfaces enable 
chemical reactions to occur that would not happen in the gas phase. 
These facts suggest that metallicity is an important parameter to 
explore in understanding star formation. However, most of what is 
known about the physics of star formation is deduced from observa- 
tions in solar-metallicity Galactic star forming regions. The nearest 
templates for the detailed study of star formation under metal-poor 
conditions are the Magellanic Clouds (MCs), with ISM metallici- 
ties of ZsMC ^ 0.2 Zq and Zlmc ^ 0.4 Z©, respectively for the 
SMC and LMC (e.g., |Russell & Dopita|1992l ). Even over the metal- 
licity range covered by the Galaxy and the MCs, cooling and chem- 
istry rates, and star formation timescales could be affected ( Banerji| 
|et al.|2009| . The low metallicity of the SMC in particular is typical 
of galaxies during the early phases of their assembly, and studies of 
star formation in the SMC provide a stepping stone to understand- 
ing star formation at high redshift where these processes cannot be 
directly observed. 

The Spitzer Space Telescope (Spitzer, [Werner et al.|2004] ) fi- 
nally allowed the identification of sizable samples of Young Stellar 
Objects (YSOs) in the MCs. The Spitzer Legacy Programs (SAGE, 
[Meixner e t aElOOG , SAGE-SMC, Gord on etal.|2011a| ) have identi- 
fied 1000s of previously unknown YSOs, both in the LMC and the 
SMC ( [Whitney et al.|2008|[GFuendl & Chu|2009l Sewilo et al in 
preparation). Follow-up spectroscopic programmes have provided 
unique insight into the abundances of ices in Magellanic YSOs, 
revealing differences in the composition of circumstellar material 
at lo wer metallicity (^Oliveira et al.|2009l [2011| [Shimonishi et al.[ 
[20T0l[Seale et al.|201 i] l 

In this paper we present a sample of 34 photometrically se- 
lected YSO candidates in the SMC, observed with Spitzer' s In- 
fraRed Spectrograph (IRS, |Houck et al.[2004| . Based on the prop- 
erties of their IRS spectra and a variety of multi- wavelength diag- 
nostics, we classify 33 sources as YSOs in the SMC (27 previously 
unknown) and 1 source as a symbiotic system. This paper is or- 
ganised as follows. After describing the sample selection and the 
different observations, the spectral properties of the SMC sources 
are analysed in detail in Section 4, with the aim of assessing the 
sources' evolutionary stage. We describe the SED properties and 
modelling, dust and ice absorption, polycyclic aromatic hydrocar- 
bon (PAH), ionic fine structure and H2 emission, and the properties 
of optical counterparts to the infrared (IR) sources. The source clas- 
sification and their evolutionary status is then discussed in Section 
5. Section 6 describes the process that allowed us to identify one of 
the sources as the first D-type symbiotic system in the SMC. 



2 DEFINING THE YSO SAMPLE 



The S^MC survey ( [Bolatto et al.|2007| ) imaged the main body of 
the SMC (central 3 square degrees) in all seven Spitzer photometric 
bands, using the Infrared Array Camera (IRAC, Fazio et al._2004j 
and the Multiband Imaging Photometer for Spitzer (MIPS, [Rieke et[ 



[al.|200 4V The resulting photometric catalogue was used to select 
the present YSO sample as described below. 

The near- and mid-IR colours of evolved stars can be as red 
as embedded YSOs (especially carbon stars, which are common in 
the SMC). This makes the 70 /xm photometry critical in identify- 
ing YSOs with cold dust (as opposed to the warm dust surround- 
ing evolved stars). Thus, we selected targets with detections in all 
IRAC bands and at 24 and 70 /xm. We also imposed a 10 mJy lower 
limit to the 8 and 24 /xm fluxes so that the sources are suitable for 
Spitzer spectroscopy. From the 160 objects with 70 /xm detections, 
we were left with 46 candidates. By cross-correlating their posi- 
tions with the SIMBAD database and the Spitzefs Reserved Ob- 
ject Catalogue, we rejected 15 objects which might be (or are near 
to) evolved stars or planetary nebulae. The 31 remaining targets 
have colours [3.6] - [5.8] > 1.35 mag and [ 8] - [24] > 2.75 mag, 
typical for embedded YSOs { [Rho et al.|2006| ). Background galaxies 
can also have colours similar to those of YSOs ( [Eisenhardt, Stern &| 
[Brodwin|2004| ). However, the selected YSO candidates are bright, 
sit well above the bulk of background galaxy contamination ( [Lee et[ 
[al.|2OO6||j0rgensen et al.|2 006 ), and most are located in active re- 
gions of star formation (Fig[TJ. Thus they are likely genuine YSOs 
and not more uniformly distributed evolved stars or background 
galaxies. YSO candidates separate better from the bulk SMC pop- 
ulation in colour-magnitude diagrams, rather than colour-colour 
diagrams ( Whit ney et al.[[2008| Sewilo et al. in preparation), as 
demonstrated in Fig.|2] Our selection criteria imply the final sample 
is composed of massive YSOs (M > lOM©). The sample is ob- 
viously not complete but it is ideally designed to identify massive 
YSOs that can be the subjects of a detailed spectroscopic analysis. 

Our strategy has proved very successful, with all but one of 
the 31 candidates confirmed as an embedded YSO or compact H II 
region (Section 4). We add to the SMC sample the three YSOs iden- 
tified by[van Loon etaLl(|2008]), MSX SMC 79, IRAS 01039-7305 
and IRAS 01042—7215, based on the appearance of their 3—4 /xm 
spectra(Section 3.3). The locations of these 34 YSO candidates 
superposed on a 70-/im MIPS image of the SMC are shown in 
Fig.[l] Fig.[2]shows [3.6] - [8.0] vs. [8.0] - [24] colour-colour and 
[8.0] vs [8.0] - [24] colour-magnitude diagrams for the 34 YSO 
candidates, together with the field SMC population ( [Gordon et aL[ 
|20TlI| ). 



3 OBSERVATIONS AND CATALOGUES 

3.1 Spitzer-IRS and MIPS-SED spectroscopy 

The mid-IR, low-resolution spectra of the 31 candidate YSOs were 
obtained using the Spitzer-IRS and the spectral energy distribution 
(SED) mode of MIPS ( [Rieke et al.|2004l ) on board Spitzer, taken 
as part of the SMC-Spec program (PID: 50240, PL: G.C. Sloan). 
This GTO programme aimed at providing a comprehensive spec- 
troscopic survey of the SMC. Scientifically its goals were to study 
dust in nearly every stage of its life cycle in the SMC in order to 
assess how the interactions of dust and its host galaxy differ from 
more metal-rich systems like the Galaxy and the LMC. The pro- 
gramme was designed to provide the SMC counterpart of the LMC 
SAGE-Spec programme ( [Kemper et al.[[2010| ), albeit at a smaller 
scale. The SMC observations targeted important classes generally 
underrepresented in the Spitzer archive of SMC sources, such as 
YSOs and compact H II regions, objects in transition to and from 
the asymptotic giant branch, and supergiants. The aforementioned 
three additional objects from [van Loon et aTjpOOS) ) were initially 
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Right Ascension 

Figure 1. SAGE-SMC MIPS 70-/im image showing the location of the 34 YSO candidates (red circles and labels); sources #01 to 31 were first identified in 
this survey (see al so Oliveira et al 12011 l[?an Loon et al. 2010b >, while sources #32 (MSX SMC 79), #33 (IRAS 01039-7305) and #34 (IRAS 01042-7215) 
were identified by |van Loon et al.| i 2008) . The star forming regions discussed in Section 4.3 are labelled in blue. 



thought to be evolved stars; their Spitzer-IRS spectra were obtained 
as part of a cycle 1 GTO programme (PID: 3277, P.L: M. Egan). 

The IRS observations of the point sources in SMC-Spec were 
performed in staring mode, using the Short-Low and Long-Low 
modules (SL and LL, respectively). The 8- and 24-/im fluxes were 
used to set exposure times, aiming at signal-to-noise ratios of ^ 60 
in SL and ^ 30 in LL. The archival spectra mentioned previously 
were also obtained in the SL and LL staring modes. All spec- 
tra were reduced following standard reduction techniques. Flat- 
fielded images from the standard Spitzer reduction pipeline (ver- 
sion SI 8. 18) were background subtracted and cleaned of "rogue" 
pixels and artefacts. Spectra were extracted individually for each 
data collection event (DCE) and co- added to produce one spectrum 
per nod position. Individual spectra are extracted both using opti- 
mal and tapered extraction; optimal extraction uses super-sampled 
point- spread function (PSF) profiles to weight the pixels from the 
spatial profiles, while tapered column extraction integrates the flux 
in a spectral window that expands with wavelength ( Lebouteiller| 
|et al.|201 C)). For each extraction, the nods were then combined to 
produce a single spectrum per order, rejecting "spikes" that appear 
in only one of the nod positions. Finally, the spectra of all the seg- 



ments were combined including the two bonus orders that are use- 
ful in correcting for discontinuities between the orders. 

While for point sources optimal extraction produces the best 
signal-to-noise ratio, tapered extraction performs better if the 
source is extended ( [Lebouteiller et al.|2010| ). Some of the sources 
in our sample could be marginally extended; we opt to use optimal 
extracted spectra but check the veracity and strength of the spectral 
features against the tapered extraction spectra. The strength of rel- 
evant emission features (Sections 4.3 and 4.4) is measured by first 
fitting a series of line segments to the continua on either side of 
feature and then integrating in Fx space. 

For 10 objects in the sample (see Table|4]) we also have MIPS- 
SED spectra, covering the wavelength range 52— 90/im, as de- 
scribed in ,van Loon et al.| ([20T0b ). The IRS and MIPS-SED Spitzer 
spectra are shown in Fig. [3] (selected examples) and Appendix B 
(complete sample). 

3.2 SAGE-SMC photometric data 

The SAGE-SMC (Surveying the Agents of Galaxy Evolution in the 
Tidally-Stripped, Low Metallicity Small Magellanic Cloud) Spitzer 
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Figure 2. Colour-colour and colour-magnitude diagrams for the YSO candidates (filled red circles), as well as the field SMC population (grey scale). 



Legacy programme (PID: 40245, RL: K. Gordon) mapped almost 
the entire SMC using IRAC and MIPS. [Gordon et al.| ( [20TTat pro- 
vide a full discussion of the observations, data reduction and cata- 
logue generation; we highlight here only the more relevant details. 

In the overlap region, the SAGE-SMC and (reprocessed) 
S^MC images were combined to produce mosaic images from 
which photometric catalogues were created (we use the "Single 
Frame -i- Mosaic Photometry" catalogue). There is a systematic off- 
set between the IRAC photometry in the SAGE-SMC and S^MC 
catalogues ( [Gordon et al.|2011a| ; since there is an excellent agree- 
ment between the SAGE-SMC fluxes and those predicted from the 
IRAC calibration stars we adopt the IRAC SAGE-SMC catalogue 
fluxes. Due to updated processing and improved calibratio ns, the 
S^MC and SAGE-SMC MIPS fluxes are also different (see Gordon 
|etal.|2011a[ for full details); we also adopt the MIPS SAGE-SMC 
catalogue fluxes. Catalogue fluxes are listed in Table [T] 

Due to the complexity of the multiple datasets, the SAGE- 
SMC catalogue relies on a stringent set of rules to create the final 
catalogues, designed to maximise reliability rather than complete- 
ness ( [Gordon et al.|2011b| ). Therefore, the SAGE-SMC catalogues 
do not provide fluxes for all the IR sources in our target list: even 
though the images clearly show a point source, a few sources are 
missing the shortest or all IRAC band fluxes for a variety of rea- 
sons. Some YSOs sit in regions where there is a complex structure 
of extended emission (characteristic of star forming environments), 
and may be also marginally extended. Other issues can arise dur- 
ing the source extraction and band merging processes causing a 
particular source not to make it into the final IRAC photometric 
catalogue, for instance variability and source confusion. As part 
of their paper on SMC YSOs identified using SAGE-SMC images 
and catalogues, Sewilo et al. (in preparation) visually inspected all 
the images of the 34 YSO candidates in this sample, and performed 
aperture photometry on the mosaic images. When SAGE-SMC cat- 
alogue fluxes are unavailable, Table[T]lists extracted aperture fluxes 
(identified by *). The adopted aperture sizes are b'' for the IRAC 
bands, and 8^^ and 15^^ respectively for MIPS 24 and 70 /xm. 



3.3 Thermal-infrared ground-based spectroscopy 

L-band (2.8—4.1 /xm) spectra of eleven bright sources in the SMC 
sample (Table [4| were obtained with the Infrared Spectrometer 
And Array Camera (ISAAC) at the European Southern Observa- 
tory (ESO) Very Large Telescope (VLT) at Paranal, on the nights of 
2006 October 28 and 29 (ESO Programme 078.C-0338, PL: J.M. 
Oliveira) to search for H2O ice absorption at 3 /xm. The standard 
IR technique of chopping and nodding was employed to cancel the 
high background. The resolving power was A/ A A ^ 500. Expo- 
sure times varied between 60 and 105 min. The hydrogen lines in 
the standard stars left remnants of at most a few per cent of the con- 
tinuum level. Telluric lines were used to calibrate the wavelength 
scale to an accuracy of A A ^ 0.002 /xm. 

M-band spectra of five bright SMC objects (Table|4]) were also 
obtained with ISAAC at the VLT, on the nights of 2009 November 
4 and 5 (ESO Programme 084.C-0955, PL: J.M. Oliveira) to search 
for the CO ice feature at 4.67 /xm. Exposure times varied between 
45 and 90 min. The M-band spectra were obtained and reduced in 
the same way as the L-band spectra. Telluric lines were used to 
calibrate the wavelength scale to an accuracy of A A ^ 0.003 /xm. 

Acquisition for the L-band spectroscopy was done using high 
spatial resolution L^-band images (A = 3.8 /xm). Magnitudes were 
obtained using aperture photometry on the targets and spectral stan- 
dard stars observed at regular intervals during each night. Magni- 
tudes were converted to fluxes using the following conversion: a 
16-mag star has a flux of 9.414 x 10~^ Jy. These resulting fluxes 
(Table [TJ were used to flux-calibrate the L-band spectra. 

Some of these spectra have been published previously; [vanj 
ILoon et al.| ( [2008) 1 first identified MSXSMC79 (source #32), 
IRAS 01039-7305 (#33) and IRAS 01042-7215 ( #34) as YSOs 
based on their L-band spectra; [Oliveira et"aL] ( |201l| ) discussed the 
ice features of sources #03, 17, 18 and 34; they did not discuss the 
spectra of source #33 since the ice detections are uncertain. 
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3.4 Near- and mid-infrared ground-based imaging 

The near-infrared imaging observations were performed with the 
SOFI (Son OF ISAAC) imager at the New Technology Telescope 
(NTT) at ESO La Silla, between 2006 October 1 and 2006 Novem- 
ber 10 in service mode (ESO Programme 078.C-0319, RL: J.Th. 
van Loon). Images were obtained with the Js and Kg filters us- 
ing the 0^^288 pixel" ^ plate scale, with a total integration time of 
1 1 min per filter. The total integration time was split into jittered 
1 min exposures to allow for efficient sky removal. For each filter, 
the jittered images were reduced and combined using standard in- 
frared reduction steps implemented within ESO's data reduction 
pipelines: detector cross-talk correction, flatfield correction, sky 
subtraction and shift- addition of jittered frames. Photometric cal- 
ibration was done using dedicated standard star observations; these 
were observed several times per night in order to estimate magni- 
tude zeropoints. Since zeropoints were well behaved, we adopted a 
single value to calibrate all observations: zj^ = 23.09 =b 0.02 mag 
and — 22.37 it 0.04 mag. PSF photometry was performed 
using the DAOPHOT package ( |Stetson|1987| ) within IRAF. Typical 
full-width at half-maximum (FWHM) of the stellar profiles were 
3.8 and 3.6 pixels, corresponding to respectively l^^l and l^^O in 
the Js- and Kg -bands. Aperture correction was performed using 
bright PSF stars. Magnitudes were converted to fluxes using the 
following conversion: a 16-mag star has fluxes of 6.112 x 10~^ Jy 
and 2.584 x 10~^ Jy, respectively in Js and Ks. Fluxes are listed in 
Table[T] In Appendix C, JsKs colour composite images are shown 
for each target (Fig.|CT]), together with SpitzerllKAC [3.6]-[5.8]-[8] 
colour composites. Some sources sit in complex cluster-like envi- 
ronments; we carefully investigated each image to identify the red- 
der source as the counterpart for the mid-IR sources (these identifi- 
cations were usually very obvious). 

SOFI imaging was obtained only for the original 31 objects, 
not for the additional objects identified by I van Loon et al.| ( |2008| ). 
However, source #32 is in the SOFI field-of-view of source #09 
and Js and Ks fluxes were also measured. For the remaining two 
sources (#33 and 34) we use JHKs fluxes and images from the IR 
Survey Facility (IRSF, |Kato"et al.|2007| ). 

The 21 brightest mid-IR objects were imaged with the VLT 
Imager and Spectrograph for the IR (VISIR) at ESO, Paranal, 
in service mode over the period 2006 September 4 to October 
3 (ESO Programme 078.C-0319, PL: J.Th. van Loon). Images 
were obtained through the narrow-band PAH2 filter (centered at 
A = 11.25 /xm, half-band width AA = 0.59 /xm). Eight standard 
stars with flux densities in the range 5 — 1 1 Jy were observed for 
photometric calibration. The plate scale was 0^^075, providing a 
field-of-view of 19" x 19" . The standard IR technique of chopping 
and nodding was employed to cancel the high background. Integra- 
tion times were 24 minutes, split in 16 exposures of 23 chop cy- 
cles each. The individual exposures were combined and corrected 
for instrumental effects using version 1.5.0 of the VISIR pipeline. 
Photometry was performed by collecting the counts within a circu- 
lar aperture centered on the zero-order maximum of the diffraction 
pattern, avoiding the first Airy ring; if the source was not detected 
upper limits were calculated. Table[T]lists the fluxes. 



servatory. More information on the instrument can be found in 
[Rodgers, Conroy & Bloxhani| ( |1988) ). Standard IRAF routines were 
used for the data reduction (bias subtraction, flat fielding, wave- 
length calibration). The blue and red spectra were joined in the in- 
terval 6160— 6300 A; the joined spectra cover the wavelength range 
3600— 9500 A. Spectral standard stars, observed with the same in- 
strumental setup as the programme objects, were used to provide 
a relative flux calibration for each object. Telluric features were 
removed by observing white dwarf standards with few intrinsic 
spectral features. The final spectra have an effective resolution of 
- 4.5 A. 

Optical spectra were obtained for 32 objects (Appendix D, 
Fig. |Dl| ). At the position of the IR source, source #24 shows no 
point source at optical wavelengths, only extended emission. There 
is clearly a point source in the Js-band (Fig. |Cl| ) but it is the faintest 
object in the sample (Table[TJ; thus it was too faint to obtain an opti- 
cal spectrum. The position of IR source #3 1 sits between two bright 
optical sources; in Fig. |Cl| it is clear that this region is very complex 
with many bright sources at short wavelengths as well as bright IR 
sources. Thus no optical spectrum was obtained for this source. 

For sources #05, 06 and 29 the optical spectra exhibit only 
Balmer absorption lines, i.e. the spectra are typical of more evolved 
stars. Since these objects are also very faint in the Js-band (Table 
[T]) it is unlikely the obtained optical spectra are associated with the 
IR source. For sources #02, 10, 11, 14, 16, 17, 18, 23 and 34 the 
spectra only show Ha emission thus we perform no further analysis 
of their optical spectra. For the remaining 20 sources we analyse the 
optical spectra in detail (Section 4.5). 

3.5.2 Radio data 

Radio free-free emission from (ultra-)compact Hll regions is a 
common signpost of massive star formation]^ Using all available 
archival data from the Australia Telescope Compact Array (ATCA) 
and the Parkes radio telescope, | Wong et al.|p011a| ) and |Crawford| 
|et al.| ( |2011| ) created new high- sensitivity, high-resolution radio- 
continuum images of the SMC at 1.42, 2.37, 4.80 and 8.64 GHz 
(wavelengths respectively 20, 13, 6 and 3 cm), as well as point 
source catalogues ( |Wong et al.|2011b]|2012| ). Of the 34 objects in 
the sample, 1 1 sources have radio-continuum detections at at least 
2 frequencies. Of these, 8 are radio point sources found within A" 
of the IR source (the positional accuracy at 20 and 13 cm). Sources 
#08 and 3 1 are in complex and extended H II regions and there- 
fore the source position is less certain; there is no point source at 
the position of source #15 but we clearly see an extended radio- 
continuum object. Radio spectral indices {Su oc v^) for these 
sources are a ^ — 0.3, consistent with their classification as com- 
pact H II regions (e.g., ,Filipovic et al.,1998^ ) — background galax- 
ies are dominated by synchrotron emission with a steeper spectral 
index. The sources detected at radio wavelengths are the most lu- 
minous in our sample (Section 4.1). 



3.5 Ancillary multi-wavelength data 

3. 5. 1 Optical spectroscopy 

Our optical spectroscopy was obtained using the Double-Beam 
Spectrograph (DBS) mounted on the Nasmyth focus of the Aus- 
tralian National University 2.3-m telescope at Siding Spring Ob- 



^ By (ultra-)compact H ii region we mean a H ii region associated with the 
early stages of massive star formation, without making any statement about 
the degree of compactness or physical properties ( ^Hoare et al.,2007j ; this is 
to distinguish these objects from "classical" evolved Hii regions. 
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3.6 Resolving YSOs at the distance of the SMC 

The Spitzer photometry has a spatial resolution of 2'\ 6'\ l^" re- 
spectively for the IRAC bands and at 24 and 70 /xm, while the res- 
olution of the IRS spectra varies between 2'/ 5 and \0" . At the dis- 
tance of the SMC (60 pc, [Szewczyk et al.|2009| ), this corresponds 
to ^ 0.6 — 5 pc for the photometry and ~ 0.8 — 3 pc for the spectra. 
For comparison the Trapezium core of the Orion Nebula cluster ex- 
tends ^ 0.5 pc. Using the near-IR photometry (FWHM typically V 
or ~ 0.3 pc) we show that the Spitzer fluxes are dominated by the 
contribution of a single bright IR source. However, it is likely that 
some IR sources encase a massive binary or a dominant massive 
star surrounded by low-mass siblings that we are unable to resolve. 
Most of our analysis focuses on the YSO chemistry and other enve- 
lope properties therefore the detailed luminosity distribution within 
the YSO source is less important. Another important issue to bear 
in mind is that the different datasets sample physically distinct spa- 
tial scales due to their different spatial resolutions (see discussion 
in Section 5.1). 



4 SPECTRAL PROPERTIES 

The spectra of YSOs are characterised by a cold dust continuum, 
and often exhibit strong silicate features at 10 and 18 /xm. Silicate 
absorption superimposed on a very red continuum is indicative of 
embedded protostellar objects (e.g., |Furlan et al.|2008| ). Ice absorp- 
tion features are another common feature in spectra of embedded 
YSOs. These objects are traditionally classified as Class I sources 
(based on their IR spectral index, |Lada|1987| ), or as Stage I sources 
(based on their modelled spectral energy distributions, |Robitaille et| 
|al.| 2006 ). As the nascent massive star becomes hotter and excites 
its environment (i.e. it develops a compact H II region), emission 
features attributed to PAHs, and atomic fine-structure and H2 emis- 
sion lines also become common. Thus the infrared spectra of YSOs 
can show a superposition of ice, dust and PAH features that can be 
difficult to disentangle, in particular in extra-galactic sources (e.g., 
[Oliveira et al.|2009| ). 

For intermediate-mass YSOs in the later stages the dust co- 
coons dissipate and become hotter, and eventually IR emission 
from a circumstellar disc dominates, with prominent silicate emis- 
sion. Such objects are usually classified as Class II or Stage II ob- 
jects. Amongst such objects are Herbig Ae/Be (HAeBe) stars, that 
are often hot enough to excite PAH emission (e.g., ,Keller et al.| 
|2008l ). 

Figures [3] and |B 1 1 show all available IR photometry and spec- 
troscopy for the 34 YSO candidates. In this section we discuss the 
main spectral features and their properties and in the next section 
we discuss object classification. 



4.1 Spectral energy distribution and dust features 

The spectral index in frequency space defined as Fu oc 
responds io /3 = —a — 1, where a is the s pectral index de- 

i log(AFA 



fined in wavelength space: a 



log A 



(Andre & Mont- 



|merle|1994] ). As determined using Spitzer photometry ot Cialactic 
YSO samples. Class I objects have a > —0.5, Class II YSOs have 
-0.5 > g > -1.8 and C lass III (disc-less YSOs) have a < -1.8 
(e.g., |Muench et al.|2007t . This corresponds io /3 < -0.5 (Class I), 
-0-5 <I3<0.8 (Class II) and finally /3 > 0.8 (Class III) sources. 
These classes are more commonly used to describe the early stages 
of intermediate- and low-mass stellar evolution. 
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Figure 3. Selected examples of SEDs of YSO candidates in the SMC. Sp 



tral indices calculated in the range 3.6—24 /im are given in each panel. The 
SEDs of all objects in the sample are provided in Appendix B. 
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We calculate the spectral indices of the objects in the sample 
between 3.6 and 24 /xm. All objects have steep spectra with /3 < 
— 1.0, consistent with a Class I classification. Closer inspection of 
the SEDs (Fig.|BTJ reveals that for #19 and 20 the SED flattens 
considerably above 20 /xm, with respectively /3 ^ —0.6 and /3 ^ 
—0.2 between 20 and 32 /xm. This suggests the presence of only 
moderate amounts of cold dust. 

As part of their survey of photometric YSO candidates, Sewilo 
et al. (in preparation) have fitted the SEDs of all the sources in 
our sample, using the grid of models and fitting tool developed by 
[Robitaille et al.|p006| 2007 ). Based on model parameters (mass 
accretion rate, stellar mass and disc mass) the fitting tool classi- 
fies objects as Stage I, II and III. The SED models do not include 
the contribution from PAH emission. Good SED fits are achieved 
for the majority of objects, and all objects are classified as Stage I 
sources. Integrated SED luminosities and stellar masses are in the 
ranges 1.5 x 10^ - 1.4 x 10^ L© and 8.2 - 30 M©; sources #20 
and 28 are the least and most luminous respectively. We should 
point out that this grid of models assumes that a single YSO source 
is responsible for the fluxes measured at different wavelengths. As 
discussed in Section 3.6, some of the SMC sources may be binaries 
or small unresolved clusters; we assume that the central source can 
be represented by the equivalent of a single luminous object. 

Prominent silicate absorption is observed in the IRS spectra of 
at least ten YSOs (indicated by / in the relevant column of Table 
|3]). A further two objects have tentative detections of silicate ab- 
sorption (/? in Table |3]): source #30 shows strong PAH emission 
while the spectrum of source #24 shows only weak features. As 
discussed by Oliveira et al.| (2009 ), without detailed modelling of 
all the spectral features in the spectrum it is often difficult to dis- 
tinguish weak silicate absorption from PAH emission complexes 
between 7 and 12 /xm. Indeed, silicate absorption is conclusively 
identified only in objects where PAH emission is relatively weak or 
non existent. The red wing of the silicate absorption band can be 
deformed by the presence of the 13-/im libration mode of H2O ice. 
This might explain the particularly extended red wing of the silicate 
feature for source #06, the source with strongest H2O and CO2 ice 
in our sample (see Section [42| . Recent modelling work ( Robinson, 
[Smith & Maldoni||2012 ) suggests that constraining envelope op- 
tical depth and dust properties is crucial in asserting whether the 
libration mode is present in a spectrum. The 18-/im silicate absorp- 
tion feature seems to be present in the spectrum of some sources, 
namely #03, 17, 18, 30, and maybe #29 as well. 

The shape of the silicate emission features is a reflection of the 
composition of the dust grains. Amorphous olivine grains create a 
typical sawtooth- shaped profile peaking at 9.8 /xm, characteristic 
of small interstellar dust grains. The presence of processed (crys- 
talline) grains is revealed by a shift in the feature's peak to longer 
wavelengths, and a significant population of larger grains broadens 
the profile (e.g., "Kessler-Silacci e t al.|2006| . However, as pointed 
out by Sargent et al. ( 2009), the presence of a population of large 
grains is difficult to identify unambiguously. 

Four objects in the sample exhibit silicates in emission (in- 
dicated by /\ in Table[3|. None of the spectra exhibits crystalline 
features. The profile of source #19 is narrow with a sharp 9.8-/im 
peak, indicating unprocessed small dust grains. This object is not 
a YSO, as detailed in Section[6] Sources #05, 14, and 20 exhibit 
broader profiles that could suggest larger grains or mixed grain pop- 
ulations. The spectra of these three objects also show PAH emission 
at 11.3 fim that, in particular for source #14, masks the underlying 
profile shape. Sources #05 and 14 exhibit a red underlying contin- 
uum while the SED of #20 is essentially flat at wavelengths longer 



than 20 /xm (see above). These objects are discussed in more detail 
in Section[5]3. 

4.2 Ice absorption bands 

In cold molecular clouds, layers of ice form on the surface of dust 
grains. H2O is by far the most abundant ice (typically 10~^ — 10~^ 
with respect to H2), followed by CO2 and CO, with a combined 
abundance of 10—30% with respect to H2O ice (e.g., van DishoeckJ 
|2004| ). Understanding ice chemistry is crucial to understanding the 
gas-phase chemistry and to probing the physical conditions within 
molecular clouds. Surface chemistry and UV and cosmic-ray pro- 
cessing of the ice mantles are thought to play an important role not 
only in the formation of more complex ice species but also H2O, 
O2 and gas-phase organic molecules. Prominent ice features in the 
IRS spectral range are found between 5 and 7 /xm (attributed to 
a mixture of H2O, NH3, CH3OH, HCOOH and H2CO ices, e.g., 
Boogert et ar]|2008|), and at 15.2 /xm (attributed to CO2 ice, e.g., 
Gerakines et al.f l999','Pontoppidan et al. 2008). The ice complex 
at 5—7 iim can be difficult to identify due to the superposition of 
numerous PAH emission features ( [Spoon et al.||2002| ). At shorter 
wavelengths, ice features of H2O and CO are found in the 3—5 /xm 
range (e.g., j01iveira et al.|201l| . Circumstellar ices have been de- 
tected in massive YSO environments in the Galaxy (e.g., [G erakines] 
let al.ll999| |Gibb et al.|2004l ) and the LMC ( |van Loon et al. 2005i 
2008, Oliv eira et al.|2009||20TT| [Scale et al.|2011||Shimonishi et| 
|al.|20101 ). 

In the SMC, H2O and CO2 ices (at 3 /xm and 15.2 /xm respec- 
tively) have been detected in the environments of sources #03, 17, 
18, 32 and 34 ( |van Loon et al. [[20081 [Oliveira et al.[[201 1). Even 
though CO ice at 4.67 /xm is detected towards YSOs in the Galaxy 
and the LMC fGibb et al. 2004J [Oliveira et al.pOl l), no CO ice 
is detected towards sources #03, 17, 18, 33 and 34 (the only SMC 
sources observed at relevant wavelengths). [Oliveira et al.[ ( |201 l[ l in- 
terpreted this as a metallicity effect. Sources #03, 06 and 17 exhibit 
a broad feature between 60—70 /im in the MIPS-SED spectrum, 
attributed to crystalline H2O ice ( [van Loon et al.|201 0b). 

In this work we newly identify ice features towards a number 
of SMC YSO candidates. We detect H2O and CO2 ice absorption 
in the environments of sources #02, 05, 06, 22 and 29 (Figs.|4]|5] 
and [6]). For sources #05 and 29 we do not have 3— 4/xm spectra 
to probe the strongest 3-/xm H2O band, but we detect the 6-/xm 
H2O band — also identified in sources #03, 06, 17 and 18 (Fig[6|. 
Additionally we detect only CO2 ice towards sources #08, 13 and 
21 (Fig.[5]); no 3-/xm spectra are available and no ice is detected at 
5—7 /xm due to the presence of strong PAH emission (Table|3]). 

The blue edge of the 3-/xm H2O ice feature is set by the Earth's 
atmospheric cut-off. Therefore we use the sources' ii^s-band mag- 
nitudes to help constrain the continuum bluewards of the feature; 
the red continuum is constrained avoiding hydrogen emission at 
3.74 /xm (Pf7) and 4.05 /xm (Bra). For the CO2 ice we use a nar- 
row wavelength interval surrounding the feature to constrain the 
continuum. Optical depth spectra are determined by subtracting a 
polynomial fitted to this local pseudo-continuum from each spec- 
trum. 

The L-band spectrum of source #06 is very red and the signal- 
to-noise ratio in the blue wing of the 3-/xm H2O feature is poor. 
Therefore we do not measure the column density directly from 
the ISAAC spectrum. Instead we scaled a low-resolution AKARI 
spectrum of a LMC YSO (SSTISAGE1CJ051449.41 -671221.5, 
^Shimonishi et al.[[2010) to match the red wing of the feature in 
#06 (as shown in Fig.[4]), and measure the column density of the 
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Figure 4. SMC sources with H2O ice (left) and CO2 ice (right) detections. Dashed (red) lines show a low-resolution AKARI spectrum of a LMC YSO 
( [Shimonishi et al.|2010j to help constrain the blue wing of the H2O ice feature. Source #30 shows H2O ice but not CO2 ice. Many sources show hydrogen 
emission at 3.74 yum (Pf7) and 4.05 /im (Bra); sources #02 and 30 exhibit PAH emission at 3.3 /xm, and #06 and 30 at 16.45 /j,m. Source #22 shows [Neiii] 
emission at 15.6 /im. 
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Figure 4 - continued The H2O ice detection for source #33 is marginal. The CO2 ice detection for source #34 is weak but significant. 
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Figure 5. Additional SMC sources with CO2 ice detections. No 3—4 fim spectrum is available for these sources. The CO2 ice profile of source #21 is distorted 
by a couple of pixels with large errors, shifting the peak position. 
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Figure 6. SMC sources with ice detections in the 5—7 jim range. This band is attributed to H2O ice at 6 jim with smaller con tributions frorn more complex 
molecules (HCOOH, H2CO and NH3), and ammonium (NHJ) with a small contribution from methanol (CH3OH) at 6.85 iim \ Boogert et al. 



2008 



This ice 

complex can only be clearly identified when the 6.22-/im PAH feature is weak or absent. Some optical depth spectra are multiplied by a factor 2 (as indicated) 
to improve feature visibility. The red line is an optical depth spectrum of a Galactic YSO for comparison ( Zasowski et al. 12009} . 
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Table 2. Column density measurements (in 10^^ molecules cm~^) 
for H2O and CO2 ices in SMC YSO candidates. For source #06 the 
H2O column density measurement (indicated by *) was performed 
by scaling a better quality spectrum (see text). For sources #05 and 
29 H2O column densities (indicated by **) are estimated from the 
6-/im feature, rather than the 3-/im feature (see text). 



source # 


N(H20) 


N(C02) 


02 


7.2±2.2 


1.6±0.2 


03 


17.7±0.7 


1.7±0.1 


05 




7.3±0.9 


06 


*46.3±7 


19.2±0.5 


08 




2.3±0.2 


13 




2.2±0.2 


17 


21.6±0.8 


2.8±0.1 


18 


22.3±1.2 


6.0±0.2 


21 




3.0±0.3 


22 


27.7±3.3 


6.0±1.0 


29 


**32±5 


6.7±0.1 


30 


4.9±1.0 


<2 


32 


18.8±1.5 


1.5±0.2 


33 


3.9±1.0 


<0.3 


34 


16.6±0.7 


1.0±0.2 



scaled spectrum. This measurement is indicated by * in Table[2] 
We validate this technique by checking that for the other sources 
the two methods provide consistent column density measurements. 
The H2O column density measurements for sources #02 and 30 are 
uncertain, since the detection is weak and PAH emission at 3.3 /xm 
fills-in the red wing of the feature. For source #33 the H2O ice de- 
tection is very tentative (4-cr detection); we do not consider it any 
further. 

The 6-/im H2O ice feature is very complex; the contributions 
from other ice species can lead to overestimates of the H2O col- 
umn density (see discussion in Oliveira et al. 2009 and references 
therein). Therefore we do not measure the H2O ice column den- 
sity directly from the 6- /am feature. Since the peak optical depth 
re correlates with the column density, we use the tq measurements 
and the column densities derived from the 3-/im feature for sources 
#03, 06, 17 and 18 to estimate the H2O column densities for sources 
#05 and 29 (assuming a simple linear correlation). These measure- 
ments are indicated by ** in Table|2] 

The red wing of the CO2 ice feature for sources #08 and 22 
is affected by weak [Nelll] emission at 15.6 /xm (see next subsec- 
tion). Weak H2O ice is detected towards source #30 but no CO2 ice 
is identified (an upper limit is listed in Table[2|. The CO2 ice de- 
tection towards source #34 is weak but significant, since clear H2O 
ice is detected. For source #33 no CO2 ice is detected. 

Calculated ice column densities are listed in Table|2] The 
adopted band strengths are 2.0 x 10^^ and 1.1 x 10^^ cm^ per 
molecule, respectively for H2O and CO2 ( |Gerakines et al.|1995| ). 
The quoted uncertainties for measurements in Table[2|do not reflect 
uncertainties in continuum determination. We have re-calculated 
CO2 column densities for sources #03, 17, 18 and 34 using op- 
timally extracted IRS spectra ( [Oliveira et al.||20l"T] used tapered- 
column extraction spectra; the optimal spectra provide better fea- 
ture contrast due to improved signal-to-noise ratio). Sources with 
stronger ice absorption tend to be the most embedded (as measured 
by the spectral index), as observed also in Galactic samples (see 
for instance [Forbrich et al.|2010| ). However, a steep spectral index 



T 1 1 1 — I — I — I 1 1 1 1 1 1 — I — I — J7 
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10 100 
N(H20) [10'' cm-'] 

Figure 7. Column densities for H2O and CO2 ices. Red and green filled 
circles represent respectively LMC and Galactic YSOs; the dashed red and 
green lines represent the estimated N(C02)/N(H20) ratios, respectively 
0.32 and 0.2 jOliveira et al.|20lT] and references therein). SMC measure- 
ments are represented by blue circles. Two fits to the SMC data, of the type 
N{C02) =mx N{R20) + A^o(H20) (m is the slope and A^o(H20) 
is the H2O column density threshold), are shown (blue dashed lines), in- 
cluding and excluding the 6- 12m H2O measurements (open blue circles). 



does not imply the presence of ice absorption. There is weak anti- 
correlation between ice column density and the source luminosity. 

[Oliveira et al.| 12011) compared the column densities for 
H2O and CO2 ices for Galactic and LMC samples. They found 
iV(C02)/iV(H20) ratios for the LMC and the Galaxy - 0.32 and 
^ 0.2 respectively, consistent with previous determinations. For 
the SMC, rather than a constant ratio, Fig.|7] suggests that there 
is a H2O column density threshold for the detection of CO2 ice, 
something not observed in either the LMC or the Galaxy. Source 
#02 is the only SMC source with N(ll20)< 1 x 10^^ cm"^ 
and a CO2 ice detection, and as explained above A^(H20) may 
be underestimated for this source. We perform linear fits to the 
SMC data, of the form N{C02) = m x N{R20) + A^o(H20) 
where m is the slope and A^o(H20) is the H2O ice column density 
threshold for the detection of CO2 ice. Excluding sources #02 (see 
above) and #30 (upper limit), the fitted slope is ^ 0.4 — 0.6 and 
A^o(H20) 1.4 - 1.55 X 10^^ cm"^ (depending on whether the 
6-/im H2O measurements are included, open blue circles in Fig.jTJ). 

In the LMC, CO2 ice column densities are enhanced with 
respect to H2O ice, while the relative CO-to-C02 abundances 
are unchanged. CO2 production could be increased due to the 
stronger UV field and/or higher dust temperatures in the LMC. 
However such harsher conditions would also destroy CO ice (the 
most volatile ice species), something that is not observed. Instead 
[Oliveira et al.[pOlT] ) suggest that H2O ice is depleted due to the 
combined effects of a lower gas-to-dust ratio and stronger UV radi- 
ation field. This would push the onset of water ice freeze-out deeper 
into the YSO envelope therefore reducing the observed column 
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density (see their Figure 3). Forming deeper in the YSO envelope, 
CO2 ice would remain unaffected. 

In the Galaxy, the -thresholds for the detection of H2O and 
CO2 ices are statistically indistinguishable (see Oliveir a etal.j2011| 
and references therein), suggesting that the two ices species are co- 
spatial in YSO envelopes. However, in the SMC the present ob- 
servations suggest a column density threshold A^o(H20)^ 1.5 x 
10^^ cm~^ for the detection of CO2 ice; even the LMC mea- 
surements are consistent with a small threshold A^o(H20)^ 3 x 
10^^ cm~^. This suggests that in metal-poor environments part of 
the envelope may have H2O ice but not CO2 ice, supporting the 
scenario proposed by Oliveira et al.| p011| . The role of reduced 
shielding in regulating the ice chemistry mentioned above is also 
supported by the non-detection (with high confidence) of CO ice 
absorption in the spectra of five SMC sources (the only sources 
observed at 4.67 /xm). 

4.3 PAH and fine structure emission 

Numerous emission features originate from the C-C and C-H 
stretching and bending modes of PAH molecules, excited by UV 
radiation. The best studied emission bands are at 6.2, 7.7, 8.6, 11. 3 
and 12.7 /xm, ubiquitous in the spectra of compact H II regions and 
planetary nebulae. Many sources also exhibit a PAH emission com- 
plex at 17 /im. Often these emission bands are accompanied by 
fine-structure emission lines such as [Siv] at 10.5 /xm, [Nell] at 
12.8 /xm, [Nelll] at 15.6 /xm, [Sill] at 34.8 /xm, and [Slll] at 18.7 
and 33.5 /xm. These lines (unresolved in the low-resolution IRS 
modes) originate from ionised gas. The presence of both PAH and 
fine-structure line emission clearly suggests the presence of an ion- 
ising source of UV radiation. 

In the sample of 34 objects, all but two exhibit PAH emis- 
sion. Even when other PAH bands are weak, the relatively isolated 
11.3-/xm band can be easily identified. Thus we use the relative 
peak strength of this feature to assess how dominant PAH emis- 
sion is in shaping the IRS spectrum. We fit a local continuum using 
spectral points to the left and right of the feature, and measure the 
peak strength of the feature with respect to the underlying contin- 
uum, Fi 1.3 /Fc. Sources with Fi 1.3 /Fc ^ 1.05 essentially show no 
PAH emission; this is the case of sources #18 and 19. For sources 
with 1.05 < F11.3/FC ^ 1.30 only the 11.3-/xm band can be eas- 
ily identified, other bands are not clearly visible (typical examples 
are sources #03 and 20). Sources with F11.3/FC > 1.3 clearly ex- 
hibit the complete zoo of PAH features (e.g., #08), and for some 
sources the spectra are completely dominated by PAH emission 
(F11.3/FC > 3, e.g., #07). The different PAH groups are Usted 
in Table[3]as X, /, /, from absent to very strong, and are used 
in the source classification in Section 5. 

Recently, |Sandstrom et al.| ( |2012) ) analysed the properties of 
PAH emission observed towards a sample of six diverse star form- 
ing regions (SFRs) in the SMC, namely N76 and N66 in the north- 
east SMC bar, three regions in the southwest bar including N22, 
and N83/84 in the SMC Wing. Our sample comprises sources in the 
same six regions, and further extends the spatial coverage by sam- 
pling for instance the region in the SMC body that connects N 66 to 
the southwest bar (see Fig.[TJ. Fig. [8] summarises PAH properties in 
the SMC. On the left we show intensity ratios for the YSO candi- 
dates (black filled circles) and the averages for the aforementioned 
SMC SFRs (red filled circles). On the right we show the intensity of 
each main band normalised to total PAH intensity. In our sample the 
1 1 .3-/xm band is strong compared to the 7.7- and 6.2-/xm bands, and 
the 8.6-/xm band is very weak compared to total PAH strength, con- 



sistent with the'Sand strom et aL] ( |201 2) results. Our analysis does 
not reveal changes to PAH emission properties that would be evi- 
dence for the YSO's irradiation of its environment. Thus the PAH 
emission observed towards at least some YSO candidates may have 
an important environmental contribution. Nevertheless, for sources 
with strong emission, total PAH intensity correlates with the lumi- 
nosity of the source. 

Model PAH band ratios (adapted from |Draine & Li|200T] as 
described by |Sandstrom et al.|2012| ) show the regions of the ratio 
diagram occupied by neutral and ionised PAHs, indicating also the 
effect of decreasing grain size. The strength of the radiation field 
also has a slight effect on PAH ratios. As proposed by ,Sandstrom| 
|et al.| ( [20T2] ), the comparison with these models suggests that SMC 
PAHs are predominantly small and neutral. 

By comparing PAH properties in the SMC and the high- 
metallicity SINGS galaxies ( [Smith et aL||2007| ), [Sandstrom et al.| 
p012^ expl ain the observed differences as a metallicity effect: they 
speculate that, in the metal-poor SMC, PAHs form preferably in 
smaller grains and mostly neutral PAHs survive in the ISM. To 
further investigate whether the observed PAH ratios are directly 
related to metallicity we extend the comparison to also include 
the measurements compiled by jGalliano et al.|f2008| for a diverse 
sample of galaxies and individual SFRs of a range of metallicities. 
These data are show in Fig.[8](left): the SINGS and Galliano sam- 
ples are represented by open blue and green circles respectively, 
and the symbol size reflects each region's gas-phase oxygen abun- 
dances, in the intervals 12 + log(0/H) = 7.7 - 8.2, 8.2 - 8.7, 
and 8.7 - 9.2. We adopt 12 + log(0/H) = 8.0 for the SMC. Both 
the SINGS and Galliano samples cover a similar metallicity range, 
with the majority of objects in each sample above 8.3 and 8.8 dex 
respectively. However, these two samples occupy different regions 
in the PAH ratio diagram (with some overlap) despite their similar 
metallicities. The SMC measurements sit in a similar region to the 
Galliano sample, despite its lower metallicity. In summary our anal- 
ysis supports the suggestion that PAHs throughout the whole SMC 
are indeed predominantly small and neutral; however the sample 
comparisons we describe do not support a simple metallicity ex- 
planation; other global environmental parameters should also play 
a role (e.g., Haynes et al.|2010| ). 

Fine- structure emission from atomic ions can be intrinsic to 
YSO sources but it can also be a result of contamination from the 
diffuse gas or nearby H II regions. Furthermore, at this resolution, 
it is difficult to separate from PAH emission, particularly for the 
12.8-/xm [Ne ll] line. While the IRS spectra of 8 objects exhibit the 
18.7- and 33.5-/xm [S III] emission lines (Table|3]), only sources #26 
and 31 show strong 8.99-/xm [Arm], 10.5-/xm [S IV] and 15.6-/xm 
[Ne III] emission lines, also present but weaker in the spectra of 
#08 and 22. Spectral contamination from the environment is likely 
given the location of many targets in the vicinity of known H II 
regions (Fig.[T]and further discussion in Section 5.1). 



4.4 H2 emission 

H2 emission is expected to be ubiquitous in YSO environments. 
However, only when the molecular gas is heated to a few hundred 
K is H2 emission observable. Both UV radiation (released by the 
accretion process, the emerging star itself, or from nearby environ- 
ment) and shocks (created as outflows interact with the quiescent 
molecular cloud) can produce warm H2 gas (for a review see Habart| 
|et al.|2005) . Extinction-corrected excitation diagrams can be used 
to diagnose the gas conditions and constrain the excitation mecha- 
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Figure 8. PAH properties towards SMC sources. Left: 11.3/7.7 intensity ratio against 6.2/7.7 ratio. The different samples are: SMC YSO candidates (this 
work, filled black circles), SMC SFRs jSandstrom et al.|2012| filled red circles), a sample of galaxies and individual SFRs ^Galliano et al.|2008, green circles) 
and finally the SINGS galaxies (Smith et al.|2007| blue circles). Symbol sizes represent three metallicity ranges: 12 + log(0/H) = 7.7 - 8.2, 8.2 - 8.7, 
and 8.7 — 9.2 dex. Solid grey lines indicate model trends for neutral and ionised PAHs of varying grain sizes adapted from |Draine & Li||200H . Right: PAH 
intensity carried in each band normalised to total PAH intensity. 



nism of massive YSOs (e.g., |van den Ancker, Tielens & Wesselius| 
|2Q00] l. 

Several emission lines due to molecular hydrogen are included 
in the IRS range. In particular pure-rotational — transitions oc- 
cur at 5.51 S(7), 6.10 S(6), 6.91 S(5), 8.03 S(4), 9.66 S(3), 12.28 
S(2), 17.03 S(l) and 28.22 /xm S(0). The S(6) and S(4) transitions 
are difficult to disentangle from the PAH emission, the S(5) line 
can be contaminated by [Arm] at 6.99 /xm, and the S/N ratio in 
the region of the S(7) line is sometimes low. Nevertheless we have 
identified H2 emission in the spectra of the majority of objects in 
our sample: for 24 sources three or more unblended emission lines 
are measured (/ in Table[3|, while for eight sources only two weak 
lines are identified (/? in Table[3]). There are two sources with no 
detectable H2 emission lines (X in Table|3]). 

Since the lines discussed here are optically thin ( [Parmar, Lacy| 
|& Achtermann|199T| ), the measured line intensities can be used to 
derive the total number of molecules and the rotational temperature. 
For excited H2 gas in local thermodynamic equilibrium (LTE) at a 
single temperature, the density A^(H2) and temperature Trot are 
described by the Boltzmann distribution (in logarithmic form): 



E(J) , ^.r^rr^ ^cB 



2 ks Trot 



(1) 



where gj, E{J) and N{J) are respectively the statistical weight, 
energy and column density of the upper Jlevel|5 B — 59.33 cm 
is the H2 molecular constant. The statistical weight is given by 
gj — 9s J + 1), where the nuclear spin weight gs = 3 for 
ort/io-states (odd J) and gs = I for para-states (even J) — we 
assume an equilibrium ortho-to-para ratio opr = 3, appropriate for 
gas with T > 300 K (e.g., |Sternberg & Neufeld|1999| ). The column 
density N(J) is derived from the measured extinction-corrected 
line intensities I (J) as follows: 

A /(J) 



N{J) =47r- 



(2) 



hcAjQ 

where Aj are the transition probabilities and Q is the beamsize in 



steradian. More details on the method and values of relevant con- 
stants and energies can be found in th e literature (e.g., | Parmar et aL] 
pW]|Bernard-Salas & Tielens |2005|[Barsony et all2"oTol L 

To estimate the beamsize Q we multiply the FWHM of the 
PSF by the slit width. Slit widths are 3'/6 and 10''6 respectively 
for SL and LL. The FWHM is measured from the two-dimensional 
PSF as a function of wavelength, ranging from l'/9 at 5.5 1 /xm to 1" 
at 28.2 fim. The PSFs are those used for the optimal extraction of 
the spectra (Section 3.1 and |Lebouteiller et al.|2 010V We choose to 
have a beamsize that increases with wavelength since some objects 
are marginally resolved, and in such complex environments at these 
large distances, the IRS spectra do sample different sized emission 
regions. 

To estimate the extinction, we use PAHFIT { [Smith et al.|2007| ) 
since it disentangles the contributions from the dust and PAH emis- 
sion; the code returns fully mixed relative extinctions that can be 
used to correct the observed line intensities. The extinction correc- 
tions are largest for the S(3) 9.66-/im line, but still usually small 
for the targets in the present sample. 

Figure|9] shows examples of H2 excitation diagrams for sev- 
eral SMC sources, open and full circles represent respectively un- 
corrected and extinction-corrected measurements. For each object 
with more than two H2 emission line measurements a straight line 
can be fitted to the data points, allowing the determination of the 
temperature and column density, using Eq. (1). It can be seen from 
the examples in Fig.[9]that the weighted column density for S(0) is 
always high compared to the measurements for other lines. This is 
to be expected since S(0) has the lowest energy, therefore it probes 
the reservoir of H2 that is too cool to excite any other transitions. 
We discuss S(0) emission later. 

Whenever possible we perform two fits; the first makes use 
of all good measurements available except S(0), the second makes 



^ In the notation - S( JO, J' indicates the lower rotational state. 



© 2012 RAS, MNRAS 000,[l]|20| 



14 Oliveira et al. 



J(0) 



T(rot) = 654±17K 
N(H^) = (3.2±0.3)x10"cm-' 
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Figure 9. H2 excitation diagrams for selected SMC YSO candidates. For sources #06 and 29, the line intensities are corrected for extinction (see text for 
details). The results of two fits are displayed when available: a fit to all measured fines (red), or a fit to only S(l) to S(3) (green). For each fit, the derived 
rotational temperature Trot and column density N{Yi2) are indicated. The dashed red line shows the contribution of a cold H2 component, with fixed 
temperature Trot = 100 K. 



use only of S(l), S(2) and S(3) (red and green lines in Fig.|9] re- 
spectively). The reason to perform the second fit to the three tran- 
sitions closest in energy level is to investigate possible deviations 
from the simple single temperature scenario, as it could be expected 
if a smaller column of hotter H2 is also present ( [van den Ancker| 
|et al.|200'0| . In general, a single-temperature Boltzmann distribu- 
tion provides good fits to the weighted column densities within 
the measurement uncertainties. Firstly this validates the adopted 
wavelength-dependent beamsizes, rather than a single beamsize. 
Furthermore the column density of the S(2) line, the only para- 
state we are able to measure, does not deviate from those of the 
ort/io-states, suggesting that the adopted LTE value of opr = 3 
is adequate. The two fits described above (see examples in Fig.[9]) 
result in temperatures that agree within l.Scr for all objects except 
for source #06 (5-cr difference). This again shows that a single tem- 
perature fit is appropriate to describe the molecular gas responsible 
for most of the emission. The exception is source #06, for which 
the excitation diagram and fitted temperatures suggest the presence 
of a hotter H2 contribution (top left panel in Fig.|9]). 

For the fits to 4 or 5 data points the median temperature 
is Trot = 594 zb 99K (range 439 - 812 K), while for the 
3-measurement fits the median temperature is Trot = 562 ± 
100 K (range 439 — 843 K). The median densities are respectively 
iV(H2) = 1.8=b0.9xl0^^cm-2(range0.6-3.8xl0^^cm-2)and 



iV(H2) = 1.6 ± 1.5 X 10^^ cm-2 (range 0.7 - 4. 8 x 10^^ cm~^). 
When compared to a sample of Galactic YSOs ( [van den Anckerj 
|1999| ) covering the same luminosity range, the rotational temper- 
atures are similar. However, the column densities for the SMC 
YSO candidates tend to be smaller by a factor three, typically 
- 1.7 X 10^^ cm"2 rather than - 5 x 10^^ cm"^ 

Neither H2 temperature or column density seem correlated 
with source luminosity for the SMC and Galactic samples. For 
sources with strong PAH emission, H2 intensity (using S(3) as a 
proxy) is largest for sources with the largest total PAH intensity. 
However, there is no distinction between the derived H2 proper- 
ties for sources with ice and silicate absorption, and PAH emission. 
Since, as discussed in the next section, these features trace the YSO 
evolution, this implies that H2 emission does not seem to correlate 
with evolutionary stage, neither in SMC nor Galactic samples (see 
also |Forbrich et al.|2010| ). 

As already mentioned, the weighted column density for S(0) 
is too large when compared to the other transitions, suggesting a 
reservoir of quiescent cold gas. Since this cold component is poorly 
constrained, we opt to fix its temperature at 100 K (e.g., |Lahuis et| 
[ar]|2010| ) and simply adjust the column density to match the ob- 
served A^(0)/^(0). We do not constrain opr, expected to be ^ 1.6 
for this temperature (Sternberg & Neufeld 1999 ). The contribution 
of this cold component, added to the main warm component, is 
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shown in Fig.|9] (dashed red hne). The fitted column densities are 
in the range A^(H2) = 0.4 - 3.8 x 10^^ cm~^ (median density 
2.3 X 10^^ cm~^). Therefore the contribution of the cold molec- 
ular gas reservoir is substantial, even though its signature is only 
observed in the S(0) emission. 

On their own the column densities and temperatures derived 
in this way cannot constrain the excitation mechanism: there is sig- 
nificant overlap in the parameter range predicted for shocked and 
photo-dissociated gas fHaba rt et al.|2005[|Bemard-Salas & Tielens| 
|2005| ), even though higher temperatures are suggestive of shocked 
gas. Other diagnostics are available: PAH emission suggests photo- 
dissociated gas while [S I] emission at 25.25 /xm indicates shocked 
gas. As already discussed, the majority of the sources in the sample 
exhibit PAH emission. We do not detect the [S I] line in any of the 
spectra. Even though it is likely that both mechanisms contribute to 
H2 excitation, the available evidence suggests that radiation is the 
dominant excitation mechanism in these SMC sources. Since H2 
emission is also excited in photo-dissociation regions (e.g.,| Hab^ 
|et al.|200 5 ), it is possible that there is an environmental contribu- 
tion to emission observed towards the YSO candidates. This may 
be true in particular for the warmer component, the cool component 
originating from the denser more shielded regions. 

4.5 Optical spectra 

In this section we discuss the optical emission-line spectra of 20 
SMC sources — for the remaining 14 sources either no optical 
spectrum could be obtained, the spectrum is not associated with 
the IR source, or only Ha emission is detected (Section 3.5.1). 
In Hll regions and their precursors surrounding young massive 
stars, common optical emission lines are due to permitted hydro- 
gen (Balmer and Paschen) and O I emission (8446 A) and numer- 
ous forbidden emission lines: [Ol] (6300 A), [Oil] (3727, 7322, 
7332 A), [Nil] (6548, 6583 A), [Sll] (6717, 6731 A), as well as 
[O III] (4363, 4959, 5007 A) and [S ill] (9068, 9530 A). If the object 
is massive enough (early 0-type source) an appreciable He^ ion- 
isation zone develops (e.g., Draine 2011) and He I recombination 
emission is detected (at 3888, 4471, 5875, 6678, 7065 A). Hydro- 
gen emission in YSOs originates both from the accretion columns 
and outflows. Forbidden line emission can also originate in the rel- 
atively low-density environments of outflowing jets or winds (e.g., 
[Wh ite & Hillenbrand 2004) but it is also observed in PDR-like en- 
vironments (Storzer & Hollenbach 2000) . Velocity information can 
distinguish between the two excitation scenarios (e.g., |Storzer &| 
|Hollenbach|2000| ) but our spectra have insufficient velocity resolu- 
tion (e.g., ~ 205 km s~^ at the position of Ha). 

Figure [dT] shows the optical spectra of the 32 SMC sources; 
the last panel shows the rich spectrum of source #26 with the emis- 
sion lines identified. We implement a classification scheme that 
relies on the detection of progressively higher excitation energy 
emission lines. Thus the classification reflects the harshness of the 
near- YSO environment. Type I objects exhibit emission from the 
Balmer and Paschen series and O I, Type II objects add collision- 
ally excited lines of [O II], [N II] and [S II]. Type III objects exhibit 
also [Olll] emission. Type IV objects add [Slll], and finally Type 
V objects show prominent He I recombination emission. The type 
breakdown of the sample of twenty sources is as follows: two Type 

I, two Type Ull (only Ha and [S II] emission, see below), four Type 

II, one Type III, two Type IV, four Type IV (a single He I line 
identified), and four Type V; the optical spectrum of source #19 is 
discussed in Section 5. 

Of the ten sources classified as Type IV or V, all sources ex- 



hibit PAH emission and seven show IR fine-structure emission. 
Most are also detected at radio wavelengths (see Table [3]). Con- 
versely, all the objects that show IR fine- structure emission, and for 
which we have an optical spectrum, are classified as Type IV— V. 
These ten sources have the highest luminosities (L ^ 10^ L©), as 
determined from SED fits (Table[3]). This builds a consistent picture 
of these objects representing more evolved YSO candidates, i.e. ul- 
tracompact H II regions. However, four of these sources do exhibit 
ice features in their IRS spectrum (see next section), suggesting 
deeply embedded objects. 

As mentioned above the optical spectra do not have enough 
resolution to investigate the origin of the emission features. How- 
ever, we have looked for broadening of the line profiles that 
would indicate infall and/or outflow activity. In fact, a number of 
sources show evidence of broadened Ha profiles: 13 sources (in- 
dicated in Table[3} have FWHM in the range 300-440 km s~\ 
Furthermore, two sources (#32 and 33) exhibit extremely broad 
profiles (FWHM> 600kms~^) with line centroids shifted to ~ 
— 200kms~^. These two sources are classified as Type I/II, since 
besides Ha only [S II] emission is detected (equally broad). The 
profiles clearly suggest an origin in optically thick winds in the en- 
vironments of these two YSO candidates. 



5 SOURCE CLASSIFICATION 

All but one of the sources in the SMC sample are YSOs, based 
on the IR spectral properties analysed in the previous section: ice 
absorption, silicate absorption or emission, PAH and H2 emission, 
red continuum and SED fitting. Boyer et al.| ( |201 1| ) tentatively pro- 
posed that nine of these sources are very dusty evolved stars, based 
on Spitzer photometric criteria. Seven of those sources show ice 
absorption and two others show silicate absorption, and are thus 
clearly SMC YSOs. The only non-YSO source in our sample (#19) 
is discussed in Section|6] 

Two YSO spectral classification schemes have been recently 
developed and applied to YSO samples in the LMC. |Seale et aT] 
( 2009) performed an automated spectral classification of YSOs us- 
ing principal component analysis to identify the spectral features 
that dominate the spectra; the following classes are relevant for the 
SMC sample: S objects show spectra dominated by silicate absorp- 
tion, P and PE objects show PAH emission and IR fine-structure 
emission and O objects show silicate emission. There is an evolu- 
tionary sequence associated with this classification in the sense that 
objects with S spectra are more embedded while P and PE spectra 
are associated with more evolved compact H II regions. We have 
classified the objects in the YSO sample using this classification 
scheme (Table[3]). 

The other classification scheme was introduced by [Woods et| 
[aL] ( |201l| ): objects are classified in groups depending on spectral 
features present in the spectrum, from Gl (ice absorption), G2 (sil- 
icate absorption), G3 (PAH emission) and G4 (silicate emission). 
The main difference from the Scale classification is the usage of ice 
absorption as a clear indicator of a very cool envelope, identifying 
the earliest embedded sources. Thus we discuss the classification 
of the SMC YSOs using this method in more detail (Table[3]). 

5.1 Embedded YSOs 

There are fourteen SMC sources classified as belonging to group 
Gl, i.e. exhibiting ice absorption. Most Gl sources (eight) also ex- 
hibit silicate absorption, but this is not always the case: source #05 
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Table 3. Detailed properties of SMC YSO candidates. We use the IRS spectra to investigate PAH, fine- structure and H2 emission and silicate absorp- 
tion/emission, and calculate the spectral index between 3.6 and 24 12m. The presence or absence of a feature is indicated by / and X, ? indicating doubt (see 
discussion in the text). For PAH emission the contrast of the I13-/j,m feature with respect to the continuum is analysed in more detail: /, / indicate very 
strong, strong and weak emission respectively (see text). Four objects with silicate emission are identified with /\. For H2 emission, /? indicates objects 
for which only two emission lines were detected, rather than three to five. Ice species investigated are H2O (at 3.1, 6 and 60 /xm), CO (4.67 /xm) and CO2 
(15.2 /im). The optical ionisation classes are defined according to the emission lines present in the spectrum: Type I objects exhibit Balmer, Paschen and O I, 
Type II objects show hydrogen, O I, [Nil], [O ii] and [S ii] emission. Type III objects show the same lines plus [O ill]. Type IV objects further add [S ill], and 
finally Type V objects show all these lines plus He I emission. Some objects exhibit a stellar absorption spectrum and others only Ho; in emission (see text for 
discussion). Objects that exhibit Ha emission broader than 300 km s"-*^ are also indicated. The next column indicates whether the object has been detected 
at radio wavelengths; * signals an extended source. The sources are classified using features in their IRS spectrum, according to two classification schemes 
previously applied to samples of LMC YSOs ( jSeale et al. [2009 [[Woods et al.|2011) . The last column provides the luminosities determined from the SED fits. 
All sources are classified as YSOs (Section 5), except for source #19 that is a D-type symbiotic system (Section 6). 
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shows a silicate emission feature (see Section 5.3), while sources 
#08, 13, 21, 22 and 30 show strong PAH emission that could mask 
weak silicate absorption — we suspect this to be the case partic- 
ularly for #30. Of the aforementioned fourteen sources, ten show 
definite H2 emission, three have weak H2 emission, and one source 
shows no H2 emission. Eight Gl sources have weak or no PAH 
emission features while six sources show strong emission (three 
of which also show fine-structure emission). Of the fourteen Gl 
sources, three have no optical spectra and seven exhibit either only 
Ha or low-excitation emission lines. Another four Gl sources, 
namely #03, 08, 13 and 22 are classified as Type IV— V based on 
the counterpart optical spectrum; except for #03, these sources also 
show strong PAH and fine- structure emission. 

As already mentioned the presence of ice and silicate absorp- 



tion features indicates the presence of an embedded YSO while 
PAH, and forbidden and fine-structure emission hint at a com- 
pact H II region. However, many sources exhibit both ice absorp- 
tion and emission features. Galactic YSOs can exhibit both PAH 
and fine structure emission and ice absorption (e.g., W3 IRS5 and 
MonR2IRS2; Gibb et al. 2004 ), and many LIVIC YSOs with CO2 
ice signatures also exhibit PAH emission ( [Oliveira et al.|2009][Seale| 
|et al.[[2011| ). The same is true for the SMC sample, with half the 
sources exhibiting mixed-property spectra. Dust and ice features 
originate from the cooler regions of the embedded YSO envelope. 
On the other hand PAH and fine-structure emission can be excited 
not only by the emerging YSO itself but also by neighbouring mas- 
sive stars, and more generally in the larger H II complexes in which 
many YSOs reside. At the LMC and SMC distances (respectively 
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- 50 and GOkpc/Ngeo w & Kanbur|2008l|Szewczyk et al.|2 009^, it 
becomes impossible to disentangle the contributions of the differ- 
ent physical environments. The spatial resolution of a IRS spectrum 
varies typically from 3'' to 10^^ corresponding to 1 to 3pc in the 
SMC, implying that very different spatial scales are sampled at dif- 
ferent wavelengths. This makes it difficult to use spectral features 
to unequivocally constrain the evolutionary stage of the object. 

Nevertheless the six Gl YSOs with ice and silicate absorp- 
tion, weak or absent PAH emission, no fine-structure emission and 
quiescent optical spectra are the more embedded, earliest YSOs 
(sources #02, 17, 18, 29, 32 and 34). Sources #28 and 33 (the two 
G2 sources) and #24 (G2/G3 source) are still relatively embedded 
(silicate is seen in absorption) but no ice absorption is detected, sug- 
gesting warmer envelopes. At these early stages the YSO has little 
influence on the physical conditions in its envelope. The remain- 
ing Gl YSOs (sources #03, 06, 08, 13, 21, 22 and 30) show both 
ice absorption and PAH emission in their spectra. These objects are 
likely more evolved since the PAH and fine-structure emission indi- 
cates that the YSO is able to emit copious amounts of UV radiation, 
while still retaining enough of its cold envelope responsible for the 
absorption features (with the possible caveat that #06, 08 and 13 
sit in particularly complex environments, see Fig.[TJ. At this stage 
a (ultra-)compact H II region is emerging. 
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Figure 10. Detail of the optical spectrum of #19, showing the prominent 
Raman emission lines. 



5.2 PAH dominated YSOs 

There are fourteen sources in the sample that we classify as G3, 
indicating that PAH emission is the noticeable feature in their steep 
IRS spectra. Five of the sources have fine-structure emission as 
well. In terms of counterpart optical spectra, there does not seem 
to be a strong relation between PAH emission and optical type 
(Section 4.5), since four sources only exhibit Ha emission, three 
sources are Type II, one source is Type III and five sources are 
Type IV— V (for one further source we have no spectrum). 

These fourteen G3 sources are the most evolved in the sample, 
with the IR spectral features clearly indicating the presence of a 
well developed compact H II region, in the process of clearing out 
the remnant dusty envelope. 



5.3 Intermediate-mass YSOs 

In this subsection we discuss the three YSOs that exhibit silicates 
in emission in their IRS spectrum, classified as G4 (sources #14 
and 20) and source #05 that also shows ice absorption. The SED 
of #20 is steep but flattens considerably above 20 /xm. In fact, 
its IRS spectrum is very similar to that of the Galactic HAeBe 
star IRAS 04101-^3103 ( [Furlan et al1[2006l ), including the weak 
PAH emission. A spectrum and photometry of the optical coun- 
terpart have been analysed by |Martayan et al.| ( |2007| identified as 
SMC5_037102); they suggest that this source is a Herbig B[e] star, 
based on the presence of numerous Fe II and [Fe ll] emission lines 
and accretion signatures in the B aimer lines. We cannot detect the 
iron lines in our low-resolution optical spectrum, but we confirm 
that Ha is broad. The fact that the continuum flux rises from the 
Js-band to 70 /xm suggests that a tenuous dusty envelope remains. 

The SEDs of sources #05 and 14 both have silicate emission 
and steep IRS spectra throughout. Source #14 shows strong PAH 
emission, while #05 shows weak PAH and H2 emission as well as 
H2O and CO2 ice absorption. |Furlan et al.| (|2008 ) describes a num- 
ber of Galactic objects that they classify as evolved Class I YSOs; 
these objects still retain a low-density envelope, tenuous enough to 



reveal the thick accretion disc and the central star. The observed 
properties of sources #05, 14 and 20 suggests the same classifica- 
tion, with #20 the more evolved of the three sources. Note that the 
evolutionary scenario discussed in the previous subsection refers to 
massive star formation, while the present discussion addresses the 
evolution of intermediate-mass stars. 



6 A D-TYPE SYMBIOTIC STAR IN THE SMC 

Source #19 is a clear oddity since its optical spectrum (FigsjlO] 
|D1| ) shows a blue continuum with broad emission lines at 6825 
and 7082 A, as well as Ha, [N ll] and [O III] emission, with weak 
He I emission at 5875 and 7065 A, and a tentative detection of He II 
emission at 4686 A. Figure 



CI 



shows that next to the IR source 
there is a bright blue star. By looking at the 2-D spectrum (before 
extraction) we see that the continuum and emission line contribu- 
tions are offset by about 2'\ in the sense that the continuum origi- 
nates from the bright blue source, while the emission lines originate 
from the IR-bright source. 

The emission features at 6825 and 7082 A are identified as Ra- 
man scattering of the O VI resonance photons at 1032 and 1038 A 
by neutral hydrogen; these lines are usually observed in the spectra 
of symbiotic binary systems (e.g., |Schmid|ri989| ). Symbiotic ob- 
jects are interacting binaries, in which an evolved giant transfers 
material to a much hotter, compact companion; according to the 
classification criteria of B elczyhski et al.| ( |2000| ), the presence of 
the Raman scattering lines and optical emission lines is enough to 
identify a symbiotic system, even if no features of the cool giant are 
found. Recently it has been proposed that massive luminous B[e] 
stars may also exhibit Raman scattering lines ( [Torres et al.|2012| ). 

Source #19 has a steep SED from 2 to 20 /xm (Fig. [3]). Above 
20 /im the IRS spectrum flattens considerably and the source be- 
comes fainter at 70 /xm. This indicates the presence of a significant 
amount of (not very cold) dust. It shows a prominent silicate emis- 
sion feature, suggestive of unprocessed dust, and no PAH emission. 
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This is consistent with this object being a dusty (D-type) symbi- 
otic star ( A ngeloni et al.||20 07 ). The optical spectrum shows no 
molecular TiO bands that would indicate the presence of the red 
giant. This is another common feature of D-type symbiotic stars: 
the reddened asymptotic giant branch star is not detected at opti- 
cal wavelengths but it reveals itself at IR wavelengths (e.g., |Cor-| 
|radi et al.^2010j . In terms of photometry, source #19 has a colour 
Js — Ks ~ 1.1 mag, and it has not been identified as variable. 
Given the presence of Raman scattering and H and He emission, its 
brightness at 70 /xm, and the properties of its IRS spectrum, we pro- 
pose that #19 is a D-type symbiotic system in the SMC. This adds 
to the six S-type (dust-less) symbiotic systems already confirmed 
in the SMC ( |Mikolajewska|2004| ). 



7 SUMMARY 

We present a multi-wavelength characterisation of 34 YSO candi- 
dates in the SMC. The target objects are bright in the lO-jnm MIPS 
band, and the selection strategy aims at excluding both evolved star 
and bright galaxy interlopers. The basis of the analysis described 
here are low-resolution IR spectra obtained with Spitzer-IRS, sup- 
ported by Spitzer photometry (IRAC and MIPS), near-IR photome- 
try, 3—5 fim spectroscopy, low-resolution optical spectroscopy and 
radio data. The objective is to confirm the YSO nature of these 
SMC sources and characterise them. We summarise here our most 
important results. 

• Of the 34 sources 33 are spectroscopically identified as YSOs 
in the SMC. This now sizable sample adds to the SMC YSOs pre- 
viously identified by | van Loon et al.| ( |2008|[2010bj ) and [Oliveira et| 
|ari( [2nTT] ). 

• One object (source #19) is identified as a D-type symbiotic 
system, based on the presence of Raman emission at 6825 and 
7082 A and nebular emission lines, as well as prominent silicate 
emission. This is the first D-type symbiotic identified in the SMC. 

• Fourteen YSOs exhibit ice absorption in their spectra. We 
analyse H2O and CO2 ice column densities; we suggest the pres- 
ence of a significant H2O column density threshold for the detec- 
tion of CO2 ice in the SMC. The observed differences between 
Galactic, LMC and SMC samples can be explained as due to metal- 
licity. 

• We analyse PAH emission, which is ubiquitous in the sam- 
ple. We confirm previous results from S andstrom et al.| ( [20T2^ who 
propose that the grains responsible for PAH emission in the SMC 
are mostly small and neutral. Based on the comparison of different 
published samples, the observed PAH properties cannot be solely 
determined by metallicity. 

• Many objects show narrow emission lines in the IRS spectra 
attributed to molecular hydrogen. Excitation diagrams constrain the 
rotational temperature and H2 column density of the bulk of the gas 
responsible to the emission. When compared to Galactic sources 
( [van den Ancker|1999| the rotational temperatures are similar, but 
the H2 column densities in the SMC are generally smaller. Photo- 
dissociation is the dominant excitation mechanism. There does not 
seem to be a clear correlation between the detection of H2 emis- 
sion and its derived properties and the evolutionary stage of the 
YSO. For most sources there is also a significant reservoir of colder 
molecular gas (Trot ^ 100 K). 

• Of the 33 YSOs in our sample, 30 cover the main stages of 
massive YSO evolution. Based on the presence of ice and sili- 
cate absorption and weak/absent PAH emission, six YSOs are still 
deeply embedded. Three other embedded sources have envelopes 



already too warm for significant amount of ice to survive on the 
dust grains. A further seven sources, while still embedded (i.e. 
showing IR molecular absorption features), already present evi- 
dence of an emerging H II region. Finally, fourteen sources exhibit 
spectra with the hallmarks of compact H II regions (strong PAH and 
fine-structure emission). Three other sources exhibit silicate emis- 
sion (one shows ice absorption as well); these sources are probably 
sources in transition between Class I and Class II, i.e. precursors to 
intermediate-mass HAeBe stars that still retain a tenuous envelope. 

Scheduled ground-based observations will target more YSOs 
with ice signatures; column density measurements and modelling 
of the ice profiles will help constrain the reason for the observed 
differences in ice properties for samples of sub-solar metallicity. 
Some of the SMC YSOs are targets of an ongoing Herschel spec- 
troscopy program that aims to constrain the role of the main cooling 
agents (gas-phase CO and water, OH, [C ll] and [O I]) and obtain 
an inventory of the species present both in gas and solid phases. 
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Table 4. Summary of spectroscopic observations available for the 
SMC sources. 
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provide evidence of the complexity of some of the YSO candidate 
environments. 



APPENDIX D: OPTICAL SPECTRA 

In this appendix we show all the optical spectra obtained as possi- 
ble counterparts of the 34 IR sources in the SMC. As described in 
Section 3.5.1, the observed optical source is not always the correct 
counterpart for the IR source. No optical counterparts were found 
for sources #05, 06, 24, 29 and 31. The spectra form the basis of 
our optical classification scheme, using several emission lines. The 
last panel of Fig. jDl [ shows the identifications of the lines discussed 
in this work for source #26. 

This paper has been typeset from a TgX/ I^T^X file prepared by the 
author. 



APPENDIX A: SUMMARY OF SPECTROSCOPIC 
OBSERVATIONS 

Table|4] lists the spectroscopic observations available for each tar- 
get. The details of the observations are described in Section 3. 



APPENDIX B: SPECTRAL ENERGY DISTRIBUTIONS 

In this Appendix we show the SEDs of the 34 sources in the SMC 
sample (Fig. |Bl| ). The photometric data are compiled in Table[T] 
and spectra were obtained at the VLT (3—5 /xm) and with Spitzer- 
IRS and MIPS-SED (Table|4}. The derived spectral indices f3 are 
indicated in each panel. 



APPENDIX C: NIR THUMBNAIL IMAGES 

In this appendix we show thumbnail images for all sources in the 
SMC sample (Fig. |Cl| ). For sources #01—32, the Jgi^s compos- 
ite is shown on the left and the Spitzer [3.6]-[5.8]-[8] composite is 
shown on the right — see Section 3 for details on these observa- 
tions. For sources #33 and 34 we show a JHK colour-composite, 
constructed using IRSF images < |Kato et al.||2007| ). These images 
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Figure Bl. SEDs of the 34 SMC sources analysed in this work, showing all available IR photometry and spectroscopy. Spectral indices calculated in the range 
3.6—24 fim are provided. 
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Figure Bl - continued 
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Figure Bl - continued 
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Figure Bl - continued 
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Figure CI. Colour-composite images for the YSO candidates in the SMC. For each object JsKs (blue/ red) and Spitzer [3.6]-[5.8]-[8] images (blue/ green/ red) 
are shown. Each image is 1^ x V across and North is to the top and East to the left. For the last 2 objects we used IRSF images and show the JHK colour- 
composite. 
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Figure CI - continued 
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Figure Dl. Optical spectra of possible optical counterparts for the IR sources in the SMC. No spectra were obtained for sources #24 and 31, and for sources 
#05, 06 and 29 the optical spectra are not associated with the IR sources. The most conspicuous spectral features are labelled in the last panel for source #26. 



© 2012 RAS, MNRAS 000,[l]|20| 



30 Oliveira et al. 




E 1.0 



6000 7000 
wavelength (A) 




6000 7000 
wavelength (A) 




6000 7000 
wavelength (A) 




6000 7000 
wavelength (A) 



Figure Dl - continued 
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Figure Dl - continued 
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Figure Dl - continued 
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Figure Dl - continued The most conspicuous spectral features are labelled for source #26. 
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